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Abstract 

Employing Monte-Carlo simulation techniques we investigate the statistical properties of equally- 
charged particles confined in an one-dimensional box trap. We detect a crossover from a crystalline 
to a cluster phase with increasing temperature. The corresponding transition temperature depends 
separately on the number of particles N and the box size L, implying non-extensivity due to 
the long-range character of the interactions. The probability density of the spacing between the 
particles exhibits at low temperatures an accumulation of discrete peaks with an overall asymmetric 
shape. Around the transition temperature it is of a Gaussian form whereas in the high temperature 
regime it obeys an exponential decay. The high temperature behaviour shows a cluster phase 
with a mean cluster size that first increases with the temperature but finally saturates. The 
crossover is clearly identifiable also in the non-linear behaviour of the heat capacity with varying 
the temperature. Possible experimental realizations are briefly addressed. 

PACS numbers: 05.20.-y,64.60.an.,64.70.-p,37.90.+j 



I. INTRODUCTION 



Within the past decade, we have witnessed enormous progress with respect to the con- 
trolled manipulation of ions. This is in particular due to the application of the dynamics 
of ionic systems to spectroscopy [U [2], the implementation of quantum simulations [31 H] 
and the realization of quantum information processors [5j 16]. In this context many trapping 
methods [7HTT] have been suggested. Among the most widely used is the Paul trap [9] 
which allows to monitor single ions [T2] and study the stability of many ion crystals [T3HTB] . 
Experimental studies of ions confined in this type of trap have shown a transition from a 
cloud of ions to a crystalline structure with decreasing temperature [16J, [TF] . 

Along with the experiments there have also been various theoretical investigations in the 
formation of Coulomb crystals in the presence of a trapping potential. An early study of the 
structure of spherical Coulomb crystals [IB] showed that particles are arranged in concentric 
spherical shells with constant inter-shell distances and a hexagonal surface structure. Fur- 
thermore, for large systems a bcc lattice is formed in the interior [19] resembling the case 
of infinite Coulomb systems (OCP) [20]. The two-dimensional case has revealed even more 
exotic phenomena. Specifically, for the cylindrically confined Coulomb lattice a structural 
phase transition with increasing linear density has been found [21] • The finite 2-D systems 
of charged particles confined in a parabolic potential or a box presents an order-disorder 
phase transition with increasing temperature [22] • The density of the particles differs for 
the two potentials being almost constant in the inner region and decreasing while moving 
outwards for the parabolic case, whereas it increases radially forming distinct shells for the 
box. 

Recently, particular interest was dedicated to the study of one-dimensional systems and 
especially cold ions confined in a corresponding harmonic potential [231428] . Analytical 
approaches have been developped using perturbation theory around the classically minimum 
energy positions [25l [26] . An interesting thermodynamic behaviour different from both the 
2-D and 3-D case has been observed due to the interplay between the long-range interactions 
and strong correlations. These observations include a deviation of thermodynamic quantities 
from extensivity, a non-uniform charge density and a structural phase transition (linear to 
zig-zag) driven by the strength of the radial potential [27] and temperature [28] . 

Following the direction of the above studies, the present paper aims at describing 
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classically the thermodynamic behaviour of equally charged particles confined to an one- 
dimensional box within a wide temperature range. Such a problem is usually treated by 
computing the partition function of the system: 



where Z KE is the part of the partition function due to the kinetic energy of the particles, 
(3 = 1/ksT, g 2 the coupling constant and 



the Coulomb potential energy. Since Z v is not analytically accessible, we will proceed here 
by utilizing Monte-Carlo simulation techniques. 

Examining the distributions of the spacing between the particles with increasing tempera- 
ture, we will observe a transition from a discrete overall asymmetric form to a continuous 
exponential one. This fact can be interpreted as an evidence for a crossover from a crys- 
talline to cluster phase at a transition temperature T c , at which this distribution acquires 
a symmetric Gaussian form. The probability of forming clusters of increasing size increases 
with the temperature and finally saturates. We will then proceed and verify the crossover 
behaviour by the temperature dependence of the heat capacity. Our results indicate that 
the transition temperature T c and consequently every thermodynamic quantity depends on 
the the number of particles N and the length L of the box in an independent manner, thus 
presenting deviations from extensivity. 

The paper is organized as follows. In Sec. II we present, in some detail, the Monte- 
Carlo methods used in our simulations. Sec. Ill contains our results and a discussion of the 
observables used to characterize the statistical mechanics of the Coulomb system. Finally, 
Sec. IV provides a brief summary of our findings commenting on a possible experimental 
realization. 

II. MONTE-CARLO APPROACHES TO THE FINITE ION CHAIN 

In order to explore the statistical properties of the finite Coulomb chain confined in 
a box we employ the Metropolis and Wang-Landau algorithms, each possessing a better 



Z(L,N,T) = Z K . E Zu, 




(1) 
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efficiency with reference to the computation of different quantities in different temperature 
regimes. In order to be self-contained, we will provide below a brief description of them, 
discussing their advantages and deficiencies and explaining how they are employed in this 
work. Before doing so, it is important to notice that for a given number N of particles in the 
box the statistical properties of the Coulomb chain are determined by a single dimensionless 
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parameter A = L ^ T . This is seen by introducing the dimensionless variables £j = ^ in 
eq. ([2]) and inserting the resulting expression in the partition function Q. In the following 
we will exclusively use the dimensionless position variables £j. As a consequence the length 
of the box can be taken without loss of generality as L — 1. 

For the Metropolis algorithm [29] we generate initially a random configuration of 
the location of the particles in the box. The configuration is sorted in an ascending order 
(0 < £1 < £2 < • • • < £/v < 1)- Then we choose for the j-th particle a new position 
< £j < 1. The efficiency of the algorithm is greater if we impose the additional restriction 

< £j < £7+1. If AE = Vc(£'i) —Vc(£i) < the new configuration is accepted, i.e. £y = 
Otherwise we accept the new configuration with a probability P given by the Boltzmann 
factor P = e _/3A£; . This procedure defines a Monte-Carlo step (MCs) and after a considerable 
number of repetitions the equilibrium state is reached. The Metropolis algorithm satisfies the 
property of detailed balance and therefore it always converges, being also as a Markov chain 
very efficient in the evaluation of the equilibrium configuration of the particles. However, 
if we have a potential with many minima the particles can be trapped in one of them and 
the required number of Monte-Carlo steps (in order to reach the global equilibrium) grows 
dramatically. Thus, usually Metropolis is not a good algorithm for very low temperatures 
where the local minima of the potential dominate the statistical properties of a system. What 
is more, a large ensemble is needed in order to obtain a good estimation of the fluctuations 
and compute quantities such as the heat capacity Cy in strongly correlated systems like the 
one considered here. 

Recently another algorithm has been proposed in order to overcome these problems, 
namely the Wang-Landau algorithm [SU]. Contrary to the canonical ensemble based 
Metropolis, this algorithm uses the concept of the microcanonical ensemble and it in- 
tends to calculate the density of states (DOS) of a system as described below. First of all we 
choose the range of accessible energies and then divide it into a number of bins M. We start, 
assuming that we have a uniform density of states g(E) i.e. g(E) = 1 for every energy bin. 
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We then proceed as follows. For the j-th particle we choose a new position < £j < 
If g(E 2 ) < g(Ei) with E 1 = Vc(£i) and E 2 = Vc(£Q we accept the new configuration. Oth- 
erwise the new configuration is accepted with a probability P = 9 g ?j^\ ■ Each time an energy 
bin is visited we update the corresponding density of states by multiplying the existing value 
by a modification factor f, i.e. g(E) = f ■ g(E). We choose / = e 1 w 2.71828... We also 
update the energy histogram H(E) = H{E) + 1. In the original version of the algorithm |30j 
for discrete systems, the above steps were repeated until a flat histogram was obtained (e.g 
min(H(E)) < 0.8 • (H(E)), with (H(E)) being the mean value of the histogram). In order 
to improve the accuracy which is of order In / we decrease the modification factor / = yff 
and repeat the procedure. When In / « 10~ 8 the density of states does not any longer evolve 
any further and the simulation is stopped. The major problem of this algorithm is that it 
does not satisfy the detailed balance condition and so its convergence can not be strictly 
proved. Furthermore the method has been suggested originally for discrete systems with 
narrow energy landscapes. Nevertheless it has been used successfully in recent calculations 
considering also continuous systems If the low energy spectrum of the system is 

complex the criterion of the flatness of the histogram is never satisfied as some energy bins 
are never visited. To overcome this problem several alternatives have been proposed [33} IM] . 
In our work we use the technique described in [31] excluding some boundary bins from the 
flatness check. The Wang-Landau algorithm has the special advantage that once the DOS 
is obtained, we can easily derive all the thermodynamic quantities (e.g. Z = J2e 9(E) el3E ) 
from this function alone. Thus, the simulation time is significantly reduced when computing 
the temperature dependence of demanding quantities like the heat capacity Cy. 

In the present study 1.5 • 10 7 MCs were required in order to reach equilibrium with the 
Metropolis algorithm and an ensemble of 200 configurations in order to obtain the mean 
energy (E) as a function of A and the number of particles N. For quantities related to the 
positions of the particles at equilibrium the size of the ensemble was increased up to 10000 
configurations to achieve better statistics. The Wang- Landau algorithm has been realized 
with 2 ■ 10 7 MCs for each value of the modification factor until In / « 10~ 8 as stated above. 
In most calculations the number of particles N = 100 is used unless stated otherwise. 

The spatially resolved density and correlation functions, as well as (E(N)) have been 
computed with the Metropolis whereas the heat capacity Cz(A) has been determined using 
the Wang-Landau algorithm. As a check, the dependence of the mean potential energy 
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(E(A)) on the parameter A has been computed with both approaches. 

III. NUMERICAL RESULTS AND DISCUSSION 

We present and describe in the following our basic results for the temperature-dependent 
behaviour and properties of the ions in the box. Let us begin with the exploration of the 
minimum energy configuration (Subsec. A) and continue with its temperature dependence 
(Subsec. B). In the subsection C which addresses the main results of this work evidence for 
the crystalline-clustering crossover and the associated phenomenology is presented. 

A. Minimum energy configuration 

Some useful remarks on the minimum energy configuration of the ions in the box can be 
obtained simply by inspection of the corresponding Coulomb potential energy. First of all 
Vc(£i,£2, • • • , 6v) IS obviously a strictly increasing function of £1 > 0) and a strictly 
decreasing function of £jv (J^ < 0). Thus, the first and the last particles always occupy 
the edges of the line segment i.e. £1 = and £jy = 1- It is obvious that the problem 
possesses a symmetry axis with respect to the centre of the box. As a result, if the number 
of particles A" is odd, the central particle is positioned at the centre of the line segment 
(£|at/2J+i = 1/2). Furthermore, for more than 3 particles the minimum energy configuration 
is not the equidistant one contrary to the 3-D case where the density of charges is uniform. 
In fact, the difference A eg £j = £f — ^f 3 with {£f } the equilibrium positions of the particles at 
zero temperature and {^ q = the equidistant positions, is a discretization of a smooth 

function, inversion symmetric w.r.t. the centre position N/2 possessing a minimal value at 
imin J$ AT/4, a maximal value at i ma x ^ 3AT/4 and being almost linear within [imin^max] 
(Fig. [I]). Due to the fixed box's length and the long-range character of the interactions, 
the particles tend to accumulate at the edges of the box leaving larger inter-space distances 
in its middle part (Fig. |2j [3] (a)). This is opposite to the case of ions confined in an 1-D 
harmonic potential J23], 26| where the length of the chain is not fixed and the particles tend to 
accumulate in the inner region. Such a behaviour could be expected as the one-dimensional 
analogue of the 2-D system presented in [22] . 
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FIG. 1. (color online). The difference of the expected equidistant position of each particle from its 
minimum energy position as a function of the index of the particle for N = 100 particles confined 
in a box with unit length L = 1. 




FIG. 2. (color online). The minimum energy configurations of particles as obtained by the 
minimization of the potential Vc(xi, X2, ■ ■ ■ , xjy) for N = 6, 10 and 20. The particles are represented 
with the (blue) full dots and the equidistant positions are marked with the (red) vertical line 
segments. 

B. Temperature dependence of the configurations and the densities of charges 

To study the temperature dependence of the observables describing the properties of 
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the charged particles in the box we have to vary the parameter A. We define j- = ksTo 
yielding r = ^ = |r as a reduced temperature. As r increases positions different from 
the minimum energy ones become accessible to the particles. At temperatures r > 10 the 
particles start to form clusters of different sizes (Fig. [3] (b) ) . Thus the form of the possible 
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particle configurations changes dramatically. However, if we consider the mean positions of 
the particles (Fig. [3] (c),(d)) over an ensemble of 10 4 configurations we observe that for the 
high temperature case (Fig. [3] (d)) the clustering is averaged out and the resulting mean 
configuration resembles very much that of the low temperature regime (Fig. [3] (c)). A slight 
difference, however, is that for high temperatures the positions seem closer to the equidistant 
ones with less accumulation at the edges of the box, which are no longer occupied. 



(a) 



(b) 



(c) 
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FIG. 3. (color online), (a-d) Particle configurations for N = 100: (a) A single configuration of 
particles at r = 10~ 3 , (b) Same for r = 10 3 , (c) The mean configuration of particles for r = 10 , 
(d) Same for r = 10 3 . 

The properties of the mean equilibrium configuration for different temperatures can be 
further explored by considering the temperature dependence of the quantity (A eg £j) and the 
density function (p(£))- As seen in Fig. |l]for temperatures r < 1 the mean positions of the 
particles are identical with the minimum energy ones (Fig. Q. Even at r = 10 2 (black line 
with circles) the form of (A e9 £i) is preserved except from a decrease in its magnitude and a 
deviation at the edges. However, for r = 10 4 the form becomes linear keeping the symmetry 
around the centre but changing by an overall sign. 

A reference to the density of the absolute positions would be of no practical use since, 
due to the finite number of particles, the problem is discrete and it would only result in a set 
of delta functions. In previous studies [231 ES, EE] there has been a particular interest in the 
quantity = — £j denoting the inter-space distance between successive particles. It 
has been shown that its inverse is a smooth function of £j and represents the density of ions 
per unit length. To improve statistics we use the ensemble average of A£j i.e (A£j). Through 
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FIG. 4. (color online). The mean difference A eg £j of the position of each particle from its expected 
equidistant position as a function of the particle's index i for temperatures r = 10~ 5 (blue circles) 
, 10~ 3 (green line), 10 -1 (magenta squares), 1 (cyan triangles), 10 2 (black line with circles) and 
10 (red stars). Note that the r = 10 _1 , 1 curves for A eg £j are on top of those for r = 10~ 5 , 10~ 3 
and therefore hardly visible. 

interpolation, we produce the density of charge (p(£)) 111 the continuum limit. Our results 
are presented in Fig. [5] for various temperatures. We clearly observe the non-uniformity at 
the outer areas (close to the edges) for temperatures r < 1 (directly opposite to that of [23J) 
and an almost uniform behaviour for high temperatures except from the region close to the 
edges where the density becomes zero (not visible in the plot). As usual, there exist large 
thermal fluctuations in the high temperature regime. 

C. Crossover from the crystalline to cluster phase 

As has been stated above the quantity A£j is very useful in the determination of the density 
of charges p(£), a continuum limit of the distribution of the equilibrium positions {£$}. In this 
subsection we present the results for the distribution of A&, p(A£j) at various temperatures 
(Fig. [6j. It turns out that this quantity contrary to p(&) is well-defined for this discrete 
system and even more, it carries significant information for its thermodynamic behaviour. 
It has been computed through the histograms of A£j for 10 4 configurations and N = 100. 
For low temperatures (Fig. [6] (a), (b)) a crystalline structure is observed as p(A£j) contains 
a clearly discrete part of peaks for small distances. Thus, the particles occupy well-defined 
positions leaving specific distances between each other. The mean value of A£j for all the 
histograms approximately coincides with the equidistant spacing (A£) ~ 0.0101 ~ jfzj- The 




FIG. 5. (color online). The density of charge as a function of the particles' positions £ for 

temperatures r = 10~ 5 (blue line) , 1 (green circles), 10 3 (black line with triangles), 10 5 (estimated 
mean red line and fluctuations red stars). The particle number is N = 100. 



most probable value (max A£j w 0.0108) is slightly larger than the mean, denoting an anti- 
clustering effect and a crystallization of ions. The peaks are broadening and merging with an 
increase of the temperature and the overall distribution becomes smoother (Fig. |6](c)). More 
positions become accessible to the particles and the most probable value decreases tending 
to the mean one. At r = 10 the form of p(A£j) becomes almost symmetric (Fig. [6] (d)) 
and its line shape resembles a Gaussian as is verified by a x 2-n t (red line). Here maxA^ rs 
0.0097 is slightly less but actually very close to (A£). In the high temperature regime we 
observe a clustering effect (Fig. [6] (d) ) right opposite to that of the low temperature case. 
The effect is strongly pronounced as can be seen by inspecting maxA^ « 0.0036 which is 
much less than the mean A£j value. This trend is enhanced for ultra-high temperatures 
(r = 10 5 ) as illustrated in Fig. [6^f). The most probable value of the spacing is almost 
zero and the distribution is well-described by an exponential (red line fit with a coefficient 
of determination R 2 = 0.9962 [35j ) . Its characteristic length scale is approximately 0.0096 
which is also close to the equidistant spacing of the particles wzj- This reflects the fact that 
for high temperatures the particles have the opportunity to occupy almost every available 
position and each particle's position is completely independent from the positions of the 
neighboring particles. The continuous change of p(A£i) implying the transition from a 
crystalline to a cluster configuration with the change of temperature is a strong indicator of 
a crossover occurring at a transition temperature r c « 10, where the distribution becomes 
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Gaussian (Fig. |6](d)). 




FIG. 6. (color online), (a-f) The distributions of A£j, p(A^i), for different temperatures: (a) 
r = 10" 5 , (b) r = 1(T 3 , (c) r = 10 -1 , (d) r = 10 (the red line is a Gaussian fit), (e) r = 10 3 , 
(f) r = 10 5 (the red line is an exponential fit). In each figure the mean value of the inter-space 
distance (A£) and its most probable value (maxA£j) are provided. 

Cluster size and temperature. We have mentioned above that for temperatures r > 10 
the particles tend to form clusters of various sizes. Let us now explore how the temperature 
affects the size of the clusters formed. In what follows we consider as cluster size 1 the size 
of a neighborhood which includes only 1 particle i.e. the case when no cluster is formed. 
Furthermore, a particle is assumed to belong to the same cluster as another one if they 
are separated by a distance smaller than the minimum inter-particle distance A£; in the 
crystalline phase. This quantity changes with the number of particles N. We observe (Fig. 
[7j) that the probability of forming a larger cluster is generally increased with the increase 
of the temperature tending to a limiting form for ultra high temperatures (dark green line 
with stars). Due to the divergence of the Coulomb potential at small inter-particle distances 
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it is not preferable for the particles, even for very high temperatures, to build clusters with 
greater probability than that dictated by the green line mentioned above. 




FIG. 7. (color online). The probability distribution of cluster sizes for different temperatures 
(semi-logarithmic scale): r = 5 (blue line with circles), r = 50 (black line with squares), r = 5 • 10 2 
(red line with triangles), r = 10 3 (cyan line with stars), r = 5- 10 3 (light green line with diamonds), 
r = 10 4 (magenta line with crosses), r = 5 • 10 4 (brown line with triangles), r = 10 5 (dark green 
line with stars). 

This behaviour is also reflected in the dependence of the mean cluster size on the tem- 
perature which is presented in Fig. [8] for the particle numbers N = 50, 100 and 200. The 
variation with the temperature is similar for the three cases being almost constant (one) for 
low temperatures, then increasing logarithmically and finally saturating at a value close to 
1.5. For N = 100 the mean cluster size starts to increase after r = 10 which, as we have 
seen, is the transition temperature of the system. For N = 50 this temperature is slightly 
less (approximately r = 4) whereas for N = 200 it is around 30. Thus, we can make the 
crude statement that the transition temperature for this finite system increases (faster than 
linearly) with the number of particles. 

Correlation functions. Concerning the correlation functions of the particles' positions 

A(m) = - «&+m» «&)> 

(where ((...)) denotes the average both over the indices i and the ensemble of configura- 
tions), we encounter (Fig. [9] (a)) a very slow decay with increasing distance m, taken as the 
magnitude of the difference of the particles' indices. The form of these functions is very 
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FIG. 8. (color online), (a-c) The mean cluster size as a function of the reduced temperature r for 
systems with different number of particles: (a) N=50, (b) N=100, (c) N=200. 

close to a quadratic polynomial. A remarkable feature of the position correlation functions 
is that temperature does not affect them contrary to the typical case for most statistical 
systems (e.g. the spin correlation functions of the classical Ising model |36j). This however, 
is in agreement with the fact that the mean configuration of the particles doesn't change 
much with the increase of temperature (Figs. [3] (c),(d)). 

On the other hand, the correlation functions of the spacing between the particles 

A{m) = «A& +m A&}} - «A& +m }> «A&}} 

exhibits a strong dependence on the temperature (Fig. [9] (b)). In particular, for low tem- 
peratures we observe a slower decay of the correlation function A(m), whereas at high 
temperatures (r > 10) the decay is in general abrupt becoming practically zero for m > 1. 

Energy dependence on N, T , L. The mean potential energy is in general a function 
of the temperature T, the length of the box L, and the number of particles N separately 
i.e. (E) = f(T,L,N). In order to study numerically this dependence we need to deal 
with dimensionless quantities. We have already introduced the reduced temperature r and 
the dimensionless positions {^}. The next step is to introduce the reduced specific mean 
energy (e) = in analogy with the dimensionless reduced temperature r, and study 

the behaviour of (e) = h(r, N) keeping one of the parameters constant and varying the 
other. Note that the T, L dependence is reduced to a dependence on the single parameter r 
which exclusively determines the behaviour of the system as discussed above. Following this 
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FIG. 9. (color online), (a-b) The normalized correlation functions A(m)/A(0) for r = 10 -5 (blue 
line), r = 1CP 1 (cyan circles), r = 10 (gray triangles with black line), r = 10 5 (red line with stars) 
and for two different quantities: (a) the absolute positions of the particles £j (b) the inter-space 
distance A£j. Note that at figure (a) all the curves for different r are on top of each other and 
therefore not distinguishable. 

procedure we have computed (e(r)) for N = 100, using both the Metropolis and the Wang- 



Landau algorithm. Obviously (Fig. 10 (a)) both algorithms give similar results, which can 
be interpreted as an evidence of convergence. However, the Wang-Landau algorithm fails to 
describe the regime of high temperatures due to the predefined limited energy space available 
for the simulation. The reduced specific mean potential energy (e) is almost constant for 
low r, as the mean energy there is fully characterized by the minimum value of Vc (^'^---£n) 
and increases almost linearly with the logarithm of the reduced temperature ((e) oc log(r)) 
for higher temperatures. Thus for low temperatures we obtain the familiar (E) ~ L^ 1 law 
for the Coulomb system. The dependence on r (on both L and T) for higher temperatures 
can be interpreted as a consequence of the dependence of the transition temperature T c on 
L as we will discuss below with the help of the heat capacity C^. 

We have also used the Metropolis algorithm to obtain the dependence of (e) on N for 
three reduced temperatures, characterized by a different qualitative behaviour, belonging to 
the low, intermediate and high temperature regime respectively: r = 10~ 2 , 10 and 10 5 . The 



diagrams of Fig. 10 (c),(d) show a linear relation on a log- log scale which indicates a power 



law relation between the involved quantities. The results of the linear fits presented below 



lead us to the following conclusions. For Fig. 10 (b) (low temperature region) the fit contains 
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FIG. 10. (color online), (a) The dependence of the dimensionless specific mean energy (e) on the 
reduced temperature r on a semi-log scale for N = 100 (Metropolis: cyan circles; Wang-Landau: 
blue line), (b-d) The dependence of the dimensionless specific mean energy (e) on N on a log-log 
scale for (b) r = 10~ 2 (circles for numerical values, blue line fit), (c) r = 10 (triangles, red line 
linear fit), (d) r = 10 5 (squares, black line linear fit). 

not only a linear but also a logarithmic term indicating a relationship (e) oc Nlog(N). This 
logarithmic correction is expected for 1-D Coulomb systems in this temperature regime [23J , 
due to the formation of a crystal and the fixed particle positions. For larger temperatures 
the relation becomes linear resulting in a characteristic exponent around 1.12 for r = 10 
and 1.04 for r = 10 5 , thus approximating 1 for high temperatures ((e) oc N). It is clear 
that since the reduced mean specific energy (e) depends on N extensivity is violated, a fact 
that is expected for systems with long-range interactions. In summary, the mean Coulomb 
energy scales roughly as (E) oc ^- which coincides with our intuition for Coulomb systems. 

Specific heat capacity cl- In order to proceed in our study of the temperature dependence 
of (e) we have also calculated the specific heat capacity cl = \&r) °^ ^ ne Coulomb gas 
confined in an one-dimensional box. Due to the fact that numerical differentiation usually en- 
counters large errors it is difficult to extract the heat capacity's values from the results of the 



Metropolis algorithm as shown in Fig. 10 (a). Thus, we have used instead the density of states 
g(e) obtained from the Wang-Landau algorithm, in order to compute the quantities (e) , (e 2 ) 
and finally the reduced specific heat capacity via the relation: |^ = -^A 2 ((e 2 ) — (e) 2 ). We 



present our results for N = 100 in Fig. [TT] using a semi-logarithmic scale. A smooth max- 
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imum occurs which is located at r « 11, thus within the temperature region in which the 
histogram of the inter-particle distance A£j takes approximately a Gaussian form (Fig. [6] 
(d)) at the transition point from the crystalline to cluster phase. This fact reinforces the 
evidence that in this temperature region there is a crossover, whose transition temperature 
T c decreases as L~ l since for fixed N = 100 we obtain r c w 10 =>- T c « ^y-. 




FIG. 11. (color online), (a) The temperature dependence of the specific heat capacity ^ for 
N = 100 on a semi-logarithmic scale. 



IV. CONCLUSIONS 



We have explored the statistical and thermodynamical properties of a system of equally 
charged particles confined in a box of length L. This finite system exhibits a crossover 
from a crystalline to a cluster phase at a transition temperature T c which depends both 
on the number of particles N and the length L of the box. Specifically, T c is a (faster 
than linear) increasing function n(N) of N, whereas it decreases with L roughly as \. The 
thermodynamic limit in our long-range interacting system is not well-defined as it strongly 
deviates from extensivity. An implementation of non-extensivity to obtain thermodynamic 
quantities and parameters that are finite is beyond the scope of the present paper but some 
steps in this direction can be found in [37]. The thermodynamic behaviour of the finite 
system is properly described by the characteristic distributions and correlations of the scaled 
relative distances A£j of the particles, whereas quantities based on the absolute positions 
of the particles are barely affected by the temperature. Concretely, the probability density 
of A£j undergoes a marked transition from a discrete accumulative to an exponential form 
as the temperature increases. In the transition regime it acquires a Gaussian form and the 

16 



corresponding heat capacity shows a maximum. For high temperatures the particles form 
clusters whose size grows with temperature and is finally saturating. 

In order to relate the results of our simulations to concrete physical setups we need to 
assign specific values to the parameters. Thus, for example, the temperature measured in 
S.I. units (T SI ) is related to the reduced temperature r through T SI = = 4 ^!j L ~ 
1.7 • 10 _5 Z 2 ^ Km, where Z is the charge of the particles. To determine the transition 
temperature for different numbers of particles one can use the dimensionless quantity G = 
^r^, which takes a fixed value G c = n ^"°^ at the critical point. Thus we find Tf 1 = 
1.7 • 10~ 4 Z 2 l ^^q ) Km. For cold atomic ions, the formation of quasi one-dimensional 
crystals with L « 1 mm [3E] leads to a transition temperature T^ 1 10~ 1 Z 2 ^ 1 ^ K. For a 
system with N w 100 and e.g Z — 1 the critical value of the temperature is T c 57 w 10" 1 fC, 
which is easily accessible experimentally. Therefore, by tuning the values of L and iV it is 
possible to study this crystalline-clustering crossover in laboratory experiments with cold 
ions assuming an 1-D trapping geometry. 

Let us now speculate about potential applications of our detected crossover and related 
phases to high energy collisional experiments of e.g. highly charged nuclei. Although the 
geometry of the high energetic colliding ions is quasi 2-D it is interesting to estimate some 
properties assuming roughly that the 1-D picture explored here is not a missleading descrip- 
tion. In nuclear systems the length scales are of order L w 10~ 15 m which for iV = 100 
nucleons leads to T^ 1 rj 10 11 K i.e. to a relative energy E c 10 MeV. The nuclear mean 
binding energy is of the order of 8MeV per nucleon resulting in a total energy of 800MeV 
for the nucleus, somewhat larger than E c . Thus the nucleons confined in the colliding nuclei 
could be in the cluster phase. Finally, we remark that the crossover physics analyzed in this 
work could also appear in one- dimensional systems showing other long-range interactions, 
such as dipolar chains [39|. 
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